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. Abstract 

a^ 

^ \ The instanton-anti-instanton pair induced asymptotics of perturbation theory 

Q ' expansion for the cross section of electron-positron pair annihilation to hadrons 

^0 ■ and hadronic width of r-lepton was found. For Nf = Nc the nonperturbative 

\Q ', instanton contribution is finite and may be calculated without phenomenological 

input. The instanton induced peturbative asymptotics was shown to be enhanced 
as {n + 10)! and in the intermediate region n < 15 may exceed the renormalon 
contribution. Unfortunately, the analysis of ~ 1/n corrections shows that for 
I n ~ 10 the obtained asymptotic expressions are at best only the order of magnitude 

■ estimate. The asymptotic series for Re+e- ^hadrons i though obtained formally for 

CD . Nf = Nc, is valid up to energies ~ 15Gev. The instanton-anti-instanton pair 

I nonperturbative contribution to Rr^hadrons blows up. On the one hand, this 

means that instantons could not be considered ab-initio at such energies. On the 
other hand, this result casts a strong doubt upon the possibility to determine the 
as from the r-lepton width. 
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1 Introduction. Instanton — Renormalon. 



The semiclassical technique for high order perturbation theory estimates was suggested 
by L.N.Lipatov [||] almost 20 years ago. The main advantage of this approach is that 
there is no need to deal with the Feynman graphs in order to find the asymptotics. Only 
the contribution of one "the most important" classical configuration in the functional 
space (instanton) is considered, which can be associated with the sum of a huge number 
of these graphs. 

Shortly after the Lipatov's paper was published another mechanism for the factorial 
growth of the series of perturbation theory was found 0, §[. It turns out that in the 
theories with the running coupling constant a very simple chain of graphs may dominate 
in the asymptotics of perturbation theory. By analogy with the instanton this chain of 
graphs was called renormalon. 
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In the present work we will calculate the instanton induced contribution to the asymp- 
totics of perturbation theory series for the cross section of electron-positron pair anni- 
hilation and hadronic width of r-lepton. As we shell see, at n < 15 the instanton 
contribution to the n-th order may be quite competitive with the renormalon contribu- 
tion. 

During the last few years the issue of the asymptotic behaviour of perturbation 
theory in QCD and QED (see e.g. has attached the renewed interest. However, 

the main attention was paid to the renormalon asymptotics. 

There are two types of renormalons - ultraviolet and infrared . The usual form of 
the ultraviolet renormalon contribution to, for example, Re+e-^hadrons is: 

Re+e-^hadrons = ^ Rn (^^) ' ~ {—b)'^n\ , (1) 

where b = — \Nf ^ 10. In QCD the series (|l|) is sign alternating and at least 
the Borel sum of the series is well defined. The problems with summation of ultraviolet 
renormalon appear in QED, where the coefficient b is negative. 

Generally speaking, the infrared renormalon contribution to the asymptotics depends 
on the process considered. For Re+e-^hadrons and Rr^hadrons the contribution of the first 
infrared renormalon to the asymptotics reads: 




All problems associated with the asymptotic series become more clear if the physical 
quantities are presented in the form of Borel integral: 

Re+e- = J exp ( M F{t)dt . (3) 

Now the singularities of the function F{t) on the complex plane are the natural sources 
of asymptotic series. The first ultraviolet renormalon in QCD (0) corresponds to the 
singularity at negative value t = —1/6. The first singularity associated with infrared 
renormalon (|^) is located twice father from the origin, but at positive value t = 2/b. 
Thus, the terms of the series (H) grow with n more slowly than those of the series (|l]). 
On the other hand, for the infrared renormalon one has to choose the correct way of the 
integration over singularity which lies on the positive real axis. Different definitions of 
the integral (j^) will lead to the different power corrections SR ~ 

The issue of infrared renormalon is directly related with the well known problem 
of Landau pole in theories with running coupling. The singularity of Borel transform 
F(t) associated with infrared renormalon coincides with Landau pole which arises 
when the exchange of soft gluon in fermion loop is taken into account (see the example 
of renormalon type graphs in any of the papers [0]-|lO|)- At this point one may say 
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that the main unsolved problems of QCD (confinement) are also hidden in the infrared 
renormalon 0. 

As we have already said, another approach for high order perturbation theory esti- 
mates based on the consideration of the specific ("classical") large fiuctuations in the 
functional space was suggested by Lipatov. The natural example of such important 
fiuctuations in QCD is the instanton-anti-instanton pair. The generic form of instanton 
induced asymptotics appears to be: 

i?„,A~(n + 4iV,)! (5) 

(here we do not show either overall numerical factor, which also may be sufficiently large, 
or the factor rC with 7 ~ 1). On the Borel plane the instanton - anti- instanton pair 
corresponds to the pole at rather distant point t = 1, but residue in the pole turns out 
to be very large ~ {Ati / asY^" ■ 



Thus, we can see that though the renormalons (|l],|D do dominate at very large n (n> 
15), the instanton-induced contribution may dominate in the intermediate asymptotics 
n = 5 -h 15. If so, the pure renormalon behavior (|l]) will never be observed in directly 
calculated terms of perturbation theory due to a strong competition with the instanton 
contribution. 

It is clear, that the exactly known 3 -h 4 terms of perturbation theory series (for 
/3-function, Re+e-^hadrons or Rr^hadrons) are much Smaller than the estimate (^) and 
the question at what number n the perturbative series could reach the full strength (^ 
is open now. The only way to answer this question is to consider the ~ 1/n correc- 
tions to the leading asymptotics. In this paper we have found (and summed up) two 
subseries of corrections to RniA- The first one behaves like (A*"^ In (n)/n)'^, the second 
- like [Nl\n.[Nc) /n)^ . The final expressions for RniA are given below (^6]), (^). The 
~ 1/n corrections tend to decrease the asymptotic predictions and thus improve the 
agreement with the "experiment". But important is not this agreement. Unfortunately, 
the corrections look like ~ Nl/n and, thus, even for n ~ 10 our asymptotic expression 
can be used at best as the estimate of the order of magnitude. 

Up to now the instanton induced asymptotics of perturbation theory for Re+e-^hadrons 



was considered only in the paper of I.I.Balitsky |[TT[]. Many of the technical methods used 
by Balitsky were useful in our work, but nevertheless, he could not find correctly either 
asymptotics of perturbation theory, or nonperturbative contribution of the instanton - 



•^In fact, the situation with ultraviolet renormalon is also not very clear. Being written more accu- 
rately than the ultraviolet renormalon contribution takes the form: 

i?„ = An^(-&)"n! , (4) 

where the constant 7 ~ 1. It has been believed for a long time that the main contribution to the 
ultraviolet renormalon comes from the graphs with exchanging of one dressed gluonic line as in the case 
of infrared renormalon. However, in the recent work [p^ Vainshtain and Zakharov showed, that the 
contribution of these graphs is small compared to the graphs with two, three etc. dressed lines. As a 
result, only the constant 7 is known now and it is absolutely unclear how to search for the overall factor 
A. 
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anti-instanton pair in the most actual case Nf = Nc- The instanton induced contribu- 
tion to the perturbative asymptotics for the simple correlator of two gluonic currents 
j ~ has also been considered by one of the authors of the present paper in 

recent work . 

Similarly to the infrared renormalon (|^) all terms of the series (^) are of the same sign. 
Nevertheless, the problem of summing the series (|^) seems not as hopeless as that of the 
renormalon. Following G.'t Hooft [Q, the author of proposed to rewrite the integral 
over the instanton-anti-instanton pair in the Borel form by considering the action as 
a collective variable. Within this approach the well-separated instanton-anti-instanton 
pair is responsible for the singular part of Borel function, while the ambiguous strongly 
interacting instanton and anti-instanton contribute to its smooth part. On the other 
hand, the best way to describe the smooth part of the Borel transform is to calculate 
exactly the few first terms of perturbative expansion. The accurate subtraction from the 
singularity of dilute gas contribution in the toy model (double well oscillator) allowed us 
to find the finite nonperturbative instanton-anti-instanton contribution |]15[. In QCD at 
Nf = Nc the Borel integral diverges only logarithmically and the total nonperturbative 
contribution from instanton-anti-instanton pair may be found by cutting the instanton 
size at p < 1/Aqcd- 

In the next section we combine some useful results concerning the instanton-antiins- 
tanton pair and the behavior of light fermions in instanton background. The correlation 
function of two electro-magnetic currents is calculated in section 3. At Nf = N^ we cal- 
culate the finite nonperturbative instanton contribution as well as the instanton induced 
asymptotics of the perturbation theory. The issue of ~ 1/n corrections to the asymp- 
totics is also considered. Unfortunately, as we have already mentioned, the asymptotic 
expression for the instanton contribution may be considered at best only as the estimate 
of the order of magnitude up to n ~ 10. On the other hand if n is not very large the 
instanton may make completely unvisible the renormalon contribution. The possibility, 
discussed in section 4, to extrapolate our asymptotic formulas for Re+e-^hadrons, found 
for Nf = Nc, up to energy ~ 10 — 15Gev may be considered as a poor consolation. 
At least at such energies the high orders of the perturbation theory are not lost at the 
background of standard power corrections. More important is the issue of instanton 
contribution to Rr^hadrons (section 5) due to numerous attempts to extract from this 
quantity the value of as{mT-). Unfortunately, here our result is negative too. Even if 
one does not take into account the asymptotics of the perturbation theory (which also 
has to be reached), the pure nonperturbative contribution to Rr-^hadrons turns out to be 
one or two orders of magnitude larger than one can admit. In spite of the fact, that at 
energy m^-c^ ~ 2Gev our result is again only the estimate of the order of magnitude, 
we have a strong doubt about the possibility to find the accurate value of from the 
r-lepton width. 
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2 The Instanton— anti-instanton pair 



The interesting physical quantities, hke, e.g., Re+e-^hadrons^ Kiay be naturally found via 
the analytical continuation of the corresponding Euclidean correlators: 

Re+e-^hadronsiq^) = UTlIm U{-q'^ - is) , (6) 

where = J^fiavors ^/^/7m^/- Therefore, we have to calculate the contribution of the 
instanton-anti-instanton pair to the correlator. As we have said above, the strongly 
interacting instanton and anti-instanton contribute only to the smooth part of the Borel 
function F (|^) and do not effect the asymptotics of perturbation theory. The singular 
part of the Borel function is saturated by the almost non-interacting pseudoparticles. 
That is why, one can obtain the reliable prediction for the asymptotics starting from 
such ill defined object as the instanton-anti-instanton pair. More concretely, the field 
configuration relevant for the large orders of perturbation theory is a small instanton 
inside of a very large anti-instanton (or vice versa). The size of small instanton is 
regulated by the internal momentum in correlator @ qpi ~ 1. The size of anti-instanton 
(as well as the distance between the centers of pseudoparticles i? ~ ^ Pi) determines, 
how close we are to the singularity on the Borel plane. 

We are interested in the instanton-anti-instanton interaction in the leading approx- 
imation. Therefore, the simple sum of instanton and anti-instanton may be used: 

A, = UAApi + UjA^Ul , (7) 

where Ua,Ui are the constant SU{Nc) matrices in the color space. By trivial gauge 
rotation one can make Uj = 1. For "small" instanton the singular gauge seems to be 
preferable: 

(X - Xi)^{{x - XjY + pj) ' ^ ^ 

where rj^^, = T°'fj^i,, and r° are the usual Pauly matrices located in the upper left 2x2 
corner (another elements vanish) of N x N matrix describing the gluonic fields. 

Before we add the anti-instanton to (^, the singularity at x = xj is pure gauge 
singularity. In order to suppress the unphysical singularities in the sum (|^), one may 
choose A^ in any regular gauge which leads to A"^{x = xj) = 0. For example, one may 
slightly rotate the BPST anti-instanton: 



V 



Near the center of instanton the matrix S[x) satisfying the above condition has the form: 

S' = e'® , Q = B^Jyx - xi)^\ C^,y{x - xi)^,{x - xi)y , (10) 

where = -j^^ = -A^{x = xj) , i?^ = {xa- xj)^ and C^,^ = C^^ is an arbitrary 
symmetric tensor. 
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After direct calculation the classical action of the instanton-anti-instanton configu- 
ration may be found with the usual dipole-dipole interaction of pseudoparticles: 



A 2 2 

= -^{l-^h} , e = ,J'f\,2 ' h = 2\TrO\' - TrOO+ , (11) 
as [R^ + Pa) 



and O is the upper left 2x2 corner of the matrix U = U^Uj (0). 

The features of light fermions are mostly sensitive to the presence of instantons. The 
Dirac operator D for each flavor of massless fermions has two eigenfunctions with 
anomalously small eigenvalues X±. As we consider the case of almost noninteracting 
pseudoparticles it is natural to search for these eigenfunction in the form of linear com- 
bination of nonperturbed zero modes of separate pseudoparticles. Explicit expressions 
for zero modes in the background of singular instanton and regular anti-instanton are: 

VI/ - 1 P' ( ^\ VI/ - ^ TT ( ^ \ 

where 0"™" = Sam/ "\/2 for a = 1, 2 and 0"™ = for a > 2, a is color index, m = 1, 2 is 
spinor index, Eam is an antisymmetric tensor, and = (=F2, t). All correlation functions 
which are of interest for us are saturated by the region |a; — x/| ~ p/. Within this region 
the anti-instanton zero mode should be modified. It is easy to verify that the spinor 
function 



vr 



.u{ ^ (13) 



is the solution of the Dirac equation Di'^ a = at x — xj ~ p/ and approaches (0) at 
\x — Xi\ > p/. 

After diagonalization of the Dirac operator IDi^a within the subspace of zero modes 
(also the identity A+ = — A_ may be useful) one gets: 

^ (pi + i?2)3/2l^^^l ^\' \TrO\ '')- ^ ' 



Because the instantons interact very slightly the nonzero modes contribution to the 
fermionic determinant factorizes. The Green function in this case also has rather simple 
form: 

S{x, y) = 5a + + - Go + 0(0 . (15) 

Here, GQ{x—y) is the bare Green function, Gj, Ga are Green functions in the background 
of separate instanton and anti-instanton and Sx is the zero mode contribution. From 
(ll^ra) one gets: 



, ^+(x)^tiy) *-(a;)^+(y) 1 f TrO+ , , TrO ] , , 
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The Green function S{x, y) may contain only the terms which convert the right fermions 
to the left and vice versa. That is why the largest, proportional to term in Sx 

vanishes. As a result, the contributions of the zero modes and of "quantum" Green 
function Gi to the correlation functions are of the same order of magnitude. The last 
two terms in (pTSl) almost cancel Ga — Gq ~ ^ in the region of interest |x — x/| ~ pi <ti pa- 
The Green function in the instanton field (for simplicity, we put xj = 0) has been found 



in H 



27r^Giix,y) 



«7m 



(XT 



TxTy 



\x\ 



p2 + {t+x){t y) _^{zT )p 



.2n 



(r+y) / 1+75 

\y\ V 2 



+ ( c.c, X ^ y) 



(17) 



where = p'^ + and z = x — y. We use the Hermitean euclidean matrices 7^ which 
satisfy {7^,7;/} = 25^,^ and, for example, the bare Green function satisfies the equation 
-ilfidf,Go = 6{x-y). 



3 The calculation of the correlation function. 



Now, at last, we can write down the expression for the correlation function 

An 

f f 



n^,= 2 |e*^^exp|^e/^}|-5:ejT,.(a;,0) + (^e;)XW5.(0)} 



X 



l^3/2^TrO\Yf^^^dxdxidxAdpidpAdU 



where 



Tf,u{x,y) 



Tr{-f^S(x,y)-f^S{y,x)} - Tr{j^Go{x,y)-f^Goiy,x)} 
Bf,{x) = Tr{-f^,S{x,x)} . 



(18) 



(19) 



The factor 2 in front of the integral in (|18D accounts for the equal contribution from small 
anti- instanton and large instanton. The factor ~ ^^'^f/'^\TrO\'^^f in the square brackets 
accounts for the contribution of almost zero modes (|T4|) to the fermion determinant. 



The instanton density reads |T8|, |T9|, ^ 

p-NcC2+NfC3 

dip) 



2tt 



2Nc 



2n 



(iVe-l)!(iV,-2)! Va.(p), 



exp 



«s(p), 



(20) 



for MS scheme Ci = 2e^/^/TT^, Ca = 1.511, C3 = 0.292, cs - C3 = 21n2 - 1/6. Up to now, 
in many papers the wrong values of C2 and C3 have been used (for which, in particular. 



C2 — C3 = 2 In 2), though the error in |18| which had been done while passing from Pauli- 



Villars to dimensional regularization was corrected in and t'Hooft in later papers 
used the correct expression for d{p) . 
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We will also use the well known two-loop formula: 



Mnf^)44.nf^)V..., (21) 



where b = fN, - and b' = fN^ - fNfN^ + Nf/N^. If Nf = N, = 3 one has 6 = 9 
and 6' = 64. 



In order to reduce the integral (|T^ to the Borel one, we have to integrate over all 
collective variables xi,xa, Pi, Pa, U, and, also, over x at fixed action (pH]), that is to say, 
at fixed value of the combination C,h. This problem may be divided into two parts. The 
integration over x and xi is rather tedious algebraical problem due to the complicated 
form of the correlation function of quark currents in the instanton field (|l^) , (|19D , (|l^ . 
But from physical point of view, the main problem is the integration over size pA and 
position R = xa — xj of the large anti-instanton. There are two competing effects here. 
First, the factor (i(p^) ~ p\ tends to make the integral over pA divergent. Second, the 
almost zero fermion modes A^^-f ~ pf^/p^A^ ( p!^ tend to suppress the contribution of 
large anti-instantons. As a result, the value of the integral ([T8| ) (as well as the validity of 
our method) depends strongly on the number of light quarks Nf. The simple dimensional 
analysis (let us note, that the integrations over dpA and (I^xa are not independent and 
Pa ^ xa — xj owing to the constraint C,h = const (0)) shows that the critical value is 
Nf = Nc- If Nf < Nc the first effect dominates and the integral (p!8D diverges at large 
Pa- Nevertheless, just in this case the well defined instanton induced asymptotics of 
perturbation theory may be extracted from (0). For calculation of the integral (|T8|) 
beyond the perturbation theory at Nf < N^ the new physical income is necessary (for 
example, one may consider the instanton liquid). The most favorable case is Nf = Nc- 
In this case, the integral over pA in (p!8D diverges only logarithmically. As a result, we are 
able not only to obtain the asymptotics of perturbation theory, but also to calculate (at 
least with the logarithmic accuracy) the finite nonperturbative instanton-anti-instanton 

pair contribution to Re+e-^hadrons and Rr^hadrons- 

If Nf > Nc, the attraction of pseudoparticles which appears owing to fermionic 
zero modes prevails. As a result, the integral is saturated by ~ -R ~ P/ and the 
approximation of almost noninteracting pseudoparticles does not work. The problem of 
instanton induced asymptotics might also have the solution in this case, but, at least, 
this solution requires the considerable modification of the method used in the present 
paper. 



With the use of (^OD , (pTD let us extract the pa dependent part from (|T 




diPA) = m/p\) - d{p) , 0(x 



l + b^ln(x) 

4:71 



(22) 



Here and below = as{q) ~ as{pi). In order to find the leading perturbative asymp- 
totics, one may assume (j){x) = 1 (as it was done in [|lll]), but the calculation of the 
correlator (|T8D beyond the perturbation theory requires the use of the function 0(x) in 
the form (^2]) . Moreover, the expansion of 0(x) in a series in (a^ log(x))'^ generates an 
important set of preasymptotic corrections ~ (log(?T,)/n)'^ to the leading asymptotics. 
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The integral over pA and R = xa — xj for a fixed value of ^ gives: 

Now we would like to present in more compact form the expressions for T^j, and -B^ 
(p!8|). Actually, the traces T^^ and B^{x)B^{ii) are of interest because the polarization 
operator 11^^ is transverse. As for B^, the situation is rather simple. Formally, the 
Green function at coincident points involved to B^ goes to infinity. Nevertheless, the 



contribution to B^ from the one-instanton Green function (17) has to be of the form: 



Tr{-i^Gi{x, x)] = {x- xi)^f{x - xi) (24) 

because [x — xj)^ is the only possible vector. As a consequence of transversness of B^: 
dfj,B^ = 0, one immediately gets f{x) = 0. Thus, the only nontrivial contribution to 
originates from zero modes: 

B,{x)B,{y) = Tr{^^S,{x,x)}Tr{^,S,{y,y)} = -^§;^iH-2D) , (25) 



TC 



47^27^2 

^ y 



\TrO\^ ' WrO\^) ' 

Here = . 

After some tedious algebraic manipulations one gets: 



y 



{y, z) (x, z)\ ( odd over r 



Ti^z^ T2 T2 \ terms 



where z^ = [x — y)^, = ^{x + y)^ . Here the coordinates x and y are measured from 
the instanton center and hence, like it was done in the papers |T2| , the integral over 
d'^xj may be replaced by the integral over d'^r 

J T,,ix, y)d^r = ^ + ^ duf^ [uup^ - (T + p-)D} , (27) 
j B,{x)B,{y)d'r = duuu^^^{H - 2D) , 

where T = p^ + uuz^ and u = 1 — u . Now it is easy to find the Fourier transforms: 

/ e^'^%^(a;, Q)d'xd'xi = 4 + V /' du \k, - ^K, ( ^] d] ,(28) 

/ e'"'B^{x)B^{Q)d^xd''xi = 2p2 { {H - 2D) , 

J Jo uu \ y'uu J 

where Ko and K2 are the McDonald functions. 
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The last integral we are to calculate is the integral over the size of the "small" 
instanton - p. As we have seen above (|^,0), the p/-dependence of the integrand (except 
for that contained in (^81) ) in the large logarithms approximation arises as the multiplier 
d?{p)/ ~ p^''"^ . From the point of view of the Operator Product Expansion the large 
instanton and anti-instanton lead to the power correction in li^y 

(g,g.-^gV) ^^^2 j d^:,Tr{Glt\x)Git\x)) . (29) 

In our formulas this power correction corresponds to the term in the first integral 
(p8|). The integral over the instanton size for such correction to Re+e-^hadrons diverges 
like ~ J dpipf~^ . However, the ~ corrections, which account for the long- wave 
vacuum fiuctuations, are natural to be considered as the part of the infrared renormalon. 
Thus, we shall omit this term. All other terms in (123) contain the exponentially decreas- 
ing McDonald functions and, hence, the p-integral converges: 

^,,-3,,.„.!!^M^^ . (30) 



K2P 



q2b-2 (26-1) 
2,-,dpdu ^ 1?^-' (6-l)![(6-2)!p(6-3)! 
uu g26-2 (26-3)! 



The coefficient 6 ~ 10 may be considered as the large parameter. Therefore, the both 
integrals ( pOD have the form / e~'^'^'' p^^dp and are saturated by p ~ 6/g . Thus, the size 
of the "small" instanton turns out to be 6 times larger than the typical wave length. In 
essence, both right hand sides of ( pOD are ~ (6/g)^^e~^^ . 

At last, collecting together (|), ([Tsl), ( p8| ) and (|30|) , we get: 



' 3 (26-3)! 

47r 



y'dt/(ie|TrOp^^e^(^^"^=)-V(0exp{^e/^}A , (31) 



where (see (p5D) 



^ = (4^ - Irl) E 4 + (2^ - 4^) (E e/)^ 



(32) 



This result coincides, up to misprints, with that obtained in [|Tl|] (of course for 0(^) = 1). 

The formula for 11 (g^) may be rewritten in the form of Borel integral if one introduces 
the new variable t = 1 — h^: 

n(g2) = Const I /t/?7|TrO|2^/A|/i|^(^-^/) X 

'e[h) r(l-t)f(^/-^^)-iexp(-— 
Jo V as } 

e(-h) j^{t - i)f exp [-"^tydt 



(33) 
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Here = 0(|1 — t\) (22). It will be recalled that formally = 1 + 0{a) , but close to 
the singularity t = 1 taking into account the function changes completely the value of 
the integrand, ^-function in (^) takes the values 1 or in accordance with the sign of 
its argument. 

As we have already said, our method does not work at Nf > Nc . In fact, if Nf > Nc 
the integral (^) over the orientations of the matrix U diverges at h = 0. The physical 
nature of the divergency is clear. We have used the Faddeev-Popov unity in order to 
fix the value of the instanton-anti-instanton interaction Smt = ^(^ — 1) (recall that our 
formulas are valid only for |t— 1| -C 1). However, if Nf > Nc the main contribution to the 
integral comes from the configurations for which although the dipole-dipole interaction 
is formally small, the pseudoparticles themselves have the comparable sizes and strongly 
overlap. Of course, one can hardly find reliable result in this situation. In particular, 
our method does not work for such an attractive problem as the calculation of the 
perturbation theory asymptotics for T Zo^hadrons (A^c = 3, Nf = 5). 

Of course, our formulas are valid if Nf < N^. However, there are just 3 light quarks 
in nature and the case Nf < Nc seems to be mostly of academic interest. Furthermore, 
for Nf < Nc the correct result for Re+e-^hadrons (up to trivial modification of single 
instanton density, see discussion after eq. (PP])) has been obtained by Balitsky |TT|. 

Therefore, we would like to consider in detail only the most interesting case Nf = Nc- 
If, in addition, Nc = 3 the formula (^) for H(g^) takes the form: 

e8/3[7!6!lV47r\i2 / ri exp(-^t) exp (-^t) \ 

Here the averaging over SU{3) group was performed numerically (see Appendix). 

The expression (|^) is enough to find the leading asymptotics of perturbative expan- 
sion for H(g) (here we put = 1): 

n(g^) = Enn(^)" , (35) 

e8/3 [716112 

n„ = - ^ / 'J 0.510(n + ll)! = -0.156(n + ll)! . 
ovr^ lb! 

We see that instanton-anti-instanton induced contribution to the perturbation theory 
series do has a huge enhancement {n + 11)!. 

Both integrals in (^) diverge at t = 1 (although the main physical problem is to 
interpret the contribution of small t to H(g^) (0)). In the configuration space these 
divergences are related to the integration over almost noninteracting instanton and anti- 
instanton. Since the divergence is only logarithmic one can try out the physical intuition 
in order to restrict the range of integration in (|3^) . The natural cut-off for is pa <C 
1/Aqcd, or, in terms of t 

/ ^ \ 2 A2 

(36) 
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If so, the nonperturbative part of (|3J) may be found explicitly. To this end let us 
supplement the first term in (0) up to the principal value integral: 

/■oo exD ( — —t) 48 47r r°° dr Att 

n(g2) ^0.510P/ dt ^ °- > (|l-t|) + — e-^ / — exp ( x) . (37) 

Jo 1 — t 85 "'(A/g)^ X as 

By the use of the logarithmic variable y = ^ ln(l/a;) in the last integral one gets (up to 
correction ~ a^) 

— 0(x) exp ( x) ^— (1 - asvr^^dy = —— . (38) 

(A/g)2 X as Jo Slas 

Note, that if one replaces (f){x) by 1, the result for the nonperturbative correction will 
be 31/9 times larger. Finally, collecting together ( ^7\ ) and (^) one finds: 

(39) 

Let us say a few words about the obtained nonperturbative correction ~ e . Effec- 
tively, the integration over t in (|5^), (|3^) may be thought as the integration over the 
size of large anti-instanton. In the logarithmic scale In(p^A) ~ as{pA)~^ and one may 
say that the integration is performed over the (inverted) running coupling constant. 
The remarkable feature of our result is that all the values of as{pA) in the whole range 
o:s{pi) < C(s{pa) <^ 1 make the comparable contribution to the nonperturbative part. 

The formula (^) correctly accounts for the high orders of the perturbation theory and 
the nonperturbative corrections. However, the first terms of the perturbative expansion 
are wrong in (pOf). The only way to eliminate this defect is to calculate explicitly the 
first 3 — 5 — 10 terms of the perturbation theory. After such modification the expression 
(^) takes the form: 







n(g2) = ^n„e.ac*(^^j - (40) 



47r 



e8/3 [7!6!]2 /4vr\i2 / .oo 4^ 4^ 

TtHtIT — 0.510 P/ dt- -exp 1 +0.0182— e-- 

671'' 16! Va,/ \ JO 1 — t a. a. 



Here, for simplicity, we put = 1 in the principal value integral (pQ]). The value of N 
should be large enough so that the perturbative series at n = will be consistent with 
the asymptotics (|35|). It is this issue that is most crucial for the practical applications 
of our method. As one can see, the asymptotic prediction (|35D turns out to be many 
times larger than the "experimentally known" first 2 — 3 terms of the perturbative series 
for /3-function or Re+e-^hadrons- Of course, sooner or later, the series will reach the 
asymptotics (PH). However, if it happens only at > 10 (and it is presumably the case). 



the issue of the comparison with the "experiment" will be only of academic interest. 
Moreover at n > 10 the renormalons (|Tp, (|^) also may become important. 

The only consistent way to find the extent of validity of the formula (^) is to consider 
the 1/n corrections to the leading asymptotics. The analogous calculation for the double 
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well oscillator has been performed in our paper [1^. It is worth mentioning that the 
accuracy of the asymptotic formula for oscillator after taking into account of the ~ 1/n 
corrections turns out to be not worse than 12% starting from n = 5. 

Even for simple quantum-mechanics the calculation of the 1/n corrections to the 
instanton induced asymptotics requires considerable efforts. The corrections of the same 
order of magnitude arose while one takes into account the corrections to the instanton- 
anti-instanton interaction (both classical ~ e~'^'^/g'^ and quantum ~ e~^) and the two- 
loop correction ~ g"^ to the single instanton density. The complete calculation of the 
similar corrections in QCD seem does not possible. Nevertheless, the corrections may 
be found which have some additional enhancement. Firstly, these are the corrections of 
the order of ~ ln{n)/n. In all formulas for the polarization operator n(g^) we studiously 
kept the function (j){x) (^) which takes into account at the two-loop level the variation 

of the coupling constant ^^(pa)- As we have shown above the use of the function 

— 

(f)\x) is necessary for the correct calculation of the nonperturbative ~ e correction 
to n(g^) (|39|). On the other hand, expanding 0(x) in the powers of as one gets the 
~ [a^ ln(l — t)\^/{l — t) corrections to the integrand (|39D . Each next term of this 
series contains the small factor but also the more strong singularity due to additional 
ln(l —t). It is easy to show that after integration over t each factor [a^ ln(Y3^)]'^ leads 



to {\Yi.{n) /n) n\. The asymptotic formula (|35[) with all the ~ (ln(n)/n) included takes 



the form: 



n„ = _0.156fl-»^V"-""- („+„).= (41) 



n 



= -0.156 (n) *^(n + ll)! . 

As we have said in the introduction, the asymptotic series (^) at very large n behave 
as n^^n^. In (^) we have shown explicitly that taking into account of the asymptotic 
freedom for the large anti-instanton (0(a;)) leads to additional factor n^. 

We wrote explicitly the factor 5A^^ in the argument of the exponent in (^) in order 
to show that, unfortunately, the correction to the leading asymptotics is of the order of 
N'^/n ~ 10/n. We may also search for the corrections of the order of ~ N'^/n enhanced 
by the additional "large" factor ln(6) = ln(3Ai"c). The correction ~ Nl\a{n)/n (^l]) 
arose after the function as{pA) was treated consistently at large p^. However, (see 
discussion after formula (|30|)), the size of the "small" instanton also is parametrically 
large compared to the characteristic wave length p/ ~ h/q. Taking into account this 
effect leads to the additional factor in (pi]): 

^ 1-^2M,>W . (42) 

Finally, we have: 



n„ = -0.156 (SATj^n (n + 11)! . (43) 



Of course, this result is valid only at n ^ A*"^. Nevertheless it may happen that n 
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Some other ~ 1/n corrections also are easy to find. For example, one may calculate 
the correction induced by the quantum part of the instanton-anti-instanton interaction. 
In fact, in the leading approximation this correction will be taken into account if one 
replaces in the dipole-dipole interaction Sint = (|TTp the running constant as{q) 

by the asi^pipA)- 

However, the complete calculation of all ~ l/n corrections is too difficult, and not 
so actual. As we see the corrections to the leading asymptotics behave like Nl/n and, 
thus, the issue of quantitative comparison of the asymptotic prediction with the exactly 
calculated terms of the perturbative expansion will be possible only at n > N"^. 



4 -Re+e ^hadrons 

Now we are able to calculate the instanton contribution to the cross section of the 
electron-positron pair annihilation into hadrons. We have already shown the well known 
formula (||) connecting Re+ e- ^hadrons and imaginary part of the polarization operator. 
The main g^-dependence of 11 comes from the factor exp(— ^) ~ g"^''. However, at 
integer h (in particular if Nf = N^) the q~'^'' term does not give rise to the imaginary 
part. The imaginary part does appear after the more gentle dependence is taken into 
account H(g^) ~ [ln(-|2)]''''?~^''- Let us note that because the imaginary part of the large 
logarithm /mln(— g^) = vr is of the order of one, the Re+e- ^hadrons at Nf = Nc acquires 
an additional power of as compared to the polarization operator H^,^. After analytical 
continuation of the expression (RPf) one gets: 



«/j7'6'l2 /47r\iV exp(-^t) 

Re^e-.hadronsiq') = ^ — ) [^■^'^^^ J, \ _^ 

636 An , 477 A 
+ T^37T^— exp( ) , (44) 



where 



28985 as as 

9 a, , , ~\ ( ^ „a 



V;(x) = (l + -^ln(x))(l + 9^1n(x))^ . (45) 



Also, it is easy to find the asymptotics of the perturbation theory together with the 
~ 1/n corrections discussed in the end of preceding section: 

Re+e-^hadrons{(l ) = ^ Rn ( ^ j ' (^^) 

i?„ = -lle*^/^ ^--^0.510 (9^n) '"(n + 10)! ^ -813 (6561n)"^ (n + 10)! . 
15! ^ ^ 

As we have already said, the "experimentally" known first terms of the perturbation 
theory series Rn are in many orders of magnitude smaller than the naive asymptotics 
{n + 4:Nc)\. As well as for the H„ (^3]) the ~ l/ra corrections to Rn (^61) tend to decrease 
the asymptotic prediction. Furthermore, if one extends the expression (|46| ) down to 
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n ~ 1, this suppression may even compensate the huge factor 10!. Let us remind, 
however, that one may beheve the formula (|i6| ) only at > 10 (or even at n ^ 10). 
Moreover, deriving the formulas (|43|j46|) we have assumed the following hierarchy of the 
~ 1/n corrections: N"^ \n.{n)/n ^ A^"^ \n{Nc)/n ^ N'^/n. But actually, as it is seen from 
(^), even at A'c = 3 the corrections ~ N'^\n.{Nf^/n are much more important than the 
corrections ~ A*"^ ln(n)/n for any reasonable n. In such situation one also may expect a 
surprise from unknown ~ / n corrections. 

Up to now we have not considered the question at what energies our asymptotic 
formulas may be used. Firstly, we calculate the functional integral by the steepest 
descent method. Hence, at least the coupling constant as{p) should be small or, in 
other words, p <^ 1/A. Moreover we consider only the case Nf = N^, and so, the 
minimal value of p is determined by the mass of c quark p^m ~ On the other hand, 
as we have shown above, the effective inverse size of the "small" instanton turns out to 
be sufficiently smaller than the external momentum q ~ b/pj. Hence, it is natural to 
expect that our asymptotic formulas are valid within the energy region: 

3Gev <t:qc< 15Gev . (47) 



Everywhere above by the coupling constant as{q) was meant the three- flavors cou- 
pling a^lq). However, at q larger than 5Gev it will be natural to express the result 
in terms of a^^q). The relation between a^^q) and 05(0') in the leading approximation 
reads: 

+l,4^Uo(l), (48) 



tt3(g) "5(9) 3 \mcmb. 

where rUc and rrih are masses of c and b quarks. As a result, instead of formula (33) we 
get 

.2^ _ /m^y/VfTryi p , exp(-gt) 



Re+e-^hadronsiq') = Const i^-^j i^ — j P dt ^ _ °^ (49) 

Here we do not account for any ~ 1/n corrections and also do not write very small at 
qc ~ 5Gev nonperturbative correction ~ 1/q^^. 



5 The calculation of the r decay width 



During the last few years a question about the possibility to extract the value of the 
strong coupling ^^(m^) from the hadronic width of the r-lepton has been actively dis- 
cussed [^-||2^. It has been shown that the "standard", obtained in the framework of 
QCD sum rules nonperturbative corrections to the r decay width, do not exceed a few 
percents |]21| . The assumption that only first terms ~ 1/^'" are important allows one to 



fix the value of the strong coupling at the m,- pole within 10% accuracy. 

The first attempt to go beyond the "standard" sum rules was made in the paper 
of M.Porrati and P.Nason |25|. However, the obtained result seems to be of only a 
methodical interest. The authors of pSf have found the single-instanton contribution 
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in an empty (perturbative) vacuum and consequently obtained a completely negligible 
result proportional to the product of the light quark masses m^jm^ms. 

Much more reasonable approach to the calculation of the instanton contribution to 



Rr^hadrons h^s bccu demonstrated in the work of I.Balitsky et al. [26|. They have tryed 



to take into account the influence on the instanton of the long-wave nonperturbative 



vacuum fluctuations. Roughly speaking, in the paper the current quark masses 
were substituted by the effective masses — > — |7r^(gg)p|. Therewith the 
instanton contribution turns out to be comparable with the "standard" nonperturbative 
contributions. However, the authors of still do not put in doubt the possibility to 
extract the value of from Rr^hadrons- 

In this section we would like to consider the contribution of the instanton-anti- 
instanton pair to the r decay width. By considering the concrete topologically trivial 
configuration, we, unlike the authors of the papers |]25| , are able to find the asymp- 
totics of the perturbation theory. As for the nonperturbative contribution, at first sight 
one may expect that our "large" anti- instanton is only one of the examples of the long- 
wave vacuum fluctuations which give rise to the quark condensate {qq) and, hence, our 



nonperturbative correction has been already included in the result of the work |^. It 
would be actually the case if our universe was built only from one or two light quarks 
{Nf < Nc). In this case the integral over the size of the large anti-instanton for the 
polarization operator (0,^) diverges and it would be natural to evaluate it via some 
phenomenological quark condensate (qq) ~ Agj^^,. However, as we have seen in two pre- 
ceding sections, in the case Nf = the integral over with the logarithmic accuracy 
comes from the whole region 1/g <^ p^i -C 1/Aqcd- In particular, it means that we have 
found explicitly the most probable long-wave background for the small instanton. As we 
will see bellow, the instanton-anti-instanton pair induced correction to Rr-^hadrons turns 
out to be much larger than the semi-phenomenological result of the paper p6|. 



The reliability of our result also requires the separate consideration. The typical size 
of "small" instanton turns out to be 6 ~ 10 times larger than 1/q. Applied to the r 
decay it means that the typical p is of the order of b/rrir ~ 1/ Aqcd to say nothing about 
the "larger" anti-instanton. In such situation we may obtain only more or less reliable 
estimate of the effect on the order of magnitude. 

Like Re+e-^hadrons, the ratio of hadronic r decay width to its leptonic width Rr^hadrons 
may be found by the analytical continuation of the correlator of the weak currents from 
the euclidean q^ region (see e.g. p5| . 



Rr-.hadrons = "6^ / ds{l - s)2[(l + 2s)U^{-sml) + H^(-sm2)] , (50) 

J|s|=l 

U,,{q') = n^{q'){q,q. - q'5,.) + H^(g')g,g. = / rfxe*^^(j+(x)j.(0)) , 

where 

= Kd«+7/.(l + 7^)c? + KsM+7/.(l + l^)s ■ (51) 

Here the first integral is taken from upper to lower edge of the cut which goes along the 
positive real axis of the complex s-plane. In order to calculate the correlator (TO), one 
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may use the formulas of the preceding sections with minimal modifications. The first 
distinction is that the weak current ( |5TD include the product of two quark operators with 
different flavors. Hence, the nonconnected graphs which in the case of the correlator of 
two electro-magnetic currents ( ]T8| ) led to the term proportional to ( e/)^, do not appear 
in ( pOD (note, in the case of electro-magnetic currents the term ~ (X^^/)^ vanishes at 
Nf = Nc = 3 due to vanishing of the sum of quark charges). The fermionic Green 
function (|16D,(|l^ anticommute with 7^. As a result, the vector- vector and axial-axial 
contributions to the correlator (^) are equal. At last, one has to make the substitution 
E/e}(= 2/3) ^ (iKdP + \Vus\^)i^ 1) while going from (|T8D to (0). Note that both 
correlators (IT^) and ( ^OD turn out to be transverse. Thus, we have n^(g^) = 3n(g^) 
(see. (§^) and n-^(g^) = 0. Now, integrating over s, one can easily flnd Rr~*hadrons in 
the leading over as approximation 

33e^/M7!6!lV47r\" / /-^ exp (-^t) 

Rr^Hadrons= ,J — O.SlOP/ dt "° V (|l-t|) (52) 

40 15! Vfts / \ "'0 1 — t 

636 Att 

H e 

28985 as J 

This expression leads to the following asymptotics of the perturbation theory : 

Rr—fhadrons 

RflT 

Just as for Re+e-^hadrons (HI), there appear certain problems while interpreting the 
instanton-anti-instanton pair induced asymptotics of the perturbation theory for Rr^hadron. 
. In fact, the asymptotics (|53|) is valid only at n ^ 10, where the renormalon contribu- 
tions (|T|),(0) are much larger than the instanton one. We may only hope that at n ~ 10 
formula (|53D does give the correct estimate of the instanton contribution on the order 
of magnitude. If so, at n ~ 10 the instantons will dominate over renormalons. In the 
most interesting case n = 3 — 5 our formulas do not work. 

In order to compare our result with the experiment, let us consider in more detail 
the pure nonperturbative ~ e~"^ term in (|5^ ) (note that this correction in terms of 
Aqcd ( PID behaves like (A/m,-)^^). The quantities i?"^ for popular values of as{mr) 
are shown in the first column of the Table. One has to compare these values with the 
experimental value Rr^hadrons = 3.56 ± 0.03 [|23| (moreover, the last quantity contains 
a large trivial part R-r-^hadrons ~ Nc and one has to extract the value ag only from the 
remainder Rr-,hadrons — 3 = 0.56 ± 0.03). As we can see, the nonperturbative correction 
turns out to be dramatically large. 

There exist, however, the procedure (used also in the work [^) which allows to 
reduce the huge discrepancy with the experiment. The regular way to improve the 
nonperturbative correction (^) is to calculate the ~ as corrections to it. Since, as will 
be shown bellow, the result is going to be deceased by 30 — 50 times, one has to sum 
an infinite series of the corrections. Undoubtedly, the exact calculation even of the first 
correction of the order of as to ( |52D is beyond our abilities. All we can do is to use the 



as_Y 
An) 

61{n + 10)\(9^n 



(53) 



35 
2n 
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dependence of the coupling constant on the instanton size which is known from 

the renorm-invariance principle. Rather weak justification for taking into account just 
these particular corrections is the fact that they are enhanced as ln(3A'"c) compared to 
the other ones. As we have said above, the typical inverse size of the "small" instanton 
is much smaller than the external momentum (p*)^^ ~ The dependence of the 

coupling constant on the instanton size in the leading one-loop approximation has been 
already taken into account when we wrote the instanton density in the form d{p) ~ p*. 
At first sight, in order to account for the two- loop effects, one has simply to add to the 
nonperturbative contribution ([5^ ) the factor (see (|42|)) : 



e-- ^ _ e-- . (54) 

as J \ In (^j / \asJ 



Unfortunately, just for r-lepton the formula (|54| ) turns out to be very unstable. Deriving 
(0), we have implied that In (j^^ ^ 1, while actually rrir ~ bA. In order to improve 

53 . 

the situation, we will try to insert the exact function [— ln(pA)]"9" in the p- integral (pO]), 



as it was done in the paper p6[. While doing so, we effectively shift the maximum of the 



integrand from the dangerous region p ~ 1/A for 53/9 ~ 6 is sufficiently large quantity. 
As the result, the main contribution to Rr^hadrons comes from the instantons of the size 
p ~ A/rriT- rather than p ~ fo/m,- . Integrating over the Borel parameter (see (P^-( pUp ) 
before the p-integration (just after the Borel integration the full multiplier [In(pA)]'^ 
is gathered), we obtain the following expression for the nonperturbative correction to 
n(g^): 




Here we have also performed the analytical continuation (|50|) — >■ — sm^ . In order to 
calculate the integral over x, one may use the identity: 




ln(p)" = lim — . (56) 



Now, with the use of (pO]), we can find the nonperturbative correction to the r decay 
width 



KZhadrons = 1™ -i CXp ( - ) X (57) 
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\ Att de J J (-s)9-^ '^^ [2uu 5"° 

After trivial integration over x, s and u one gets: 

9e^/^ ( A-K An 

KZhadrons = „ ^1 TTZTT GXp ( - , ) (58) 



e^o TT ■ 70 ■ 31 \as{mT-) J ^^("^r) 
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Att dej ' ' 22er(18 - 2e) 

This expression also may be presented as a series in powers of a, . Note, that the 
first term of the series, corresponding to (d/de)^, is vanishing. This leads to loss of 
one factor l/as{mr) in Rr^hadrons (similarly, as in the case of Re+e-^hadrons) compared 
to the correlation function. Coefficients of the expansion of ( [58| ) in powers of a<j were 
found numerically. The resulting quantities R*JZhadrons shown in the last column 

,53 

of the Table. The terms of series in a,, generated by f 1 + , are maximal at 

n ~ 2 3 and become negligible at n > 8 . As we can see, taking into account of the 
two- loop dependence of as{p) allows to reduce the nonperturbative correction almost by 
two orders of magnitude. Nevertheless, the formula (BBf) turns out to be surprisingly 



stable. We have used also the Stirling formula for approximate calculation of the F- 
functions in (^81) and the approximate values of R*JZhadrons found in this way almost 



*hadron. 

coincide with results of numerical calculations listed in the table. 

It is seen from the table that passing from a, = 0.28 to a, = 0.29 our correction 
changes the sign. Moreover, the ratio R* / R in the interval 0.15 < a, < 0.43 changes the 
sign three times, being everywhere smaller than 1/30. Compared to the simple formula 
(0) such rich behaviour provides us one more evidence that effectively we work at very 
low energies and listed in the table result is, at best, only the estimate on the order of 
magnitude. 



Table 



asirrtr) 


^T^hadrons 


^T^hadrons 


0.28 


5.65 


-0.0229 


0.29 


17.44 


0.0507 


0.30 


49.23 


0.476 


0.32 


311.1 


6.70 


0.34 


1514.1 


44.79 


0.36 


5943.4 


190.1 



The full nonperturbative contribution of the instanton-anti-instanton pair to the hadronic r 
decay width R*^hadrons ^"^^ contribution accounting for only the leading over a, term in 
the sum (H) li^'^hadrons ^* various values of as{mr). 



One should remember that the obtained results should be compared with the exper- 
imental value {Rr^hadrons — 3)exp = 0.56 ± 0.03 . Wc scc that, in spite of all our effort, 
even for a, = 0.29 the instanton contribution to Rr-^hadrons still remains large. In this 
situation the only way out may be to ignore completely the instanton contribution to 
Rr^hadrons ■ For example, one may say that the series of the power corrections ~ 1/m" 
is also asymptotic and it is natural to cut off it somewhere at n ~ 4 — 8 (the instanton 
contribution behaves like ~ 1/m}^). Nevertheless, we do not know to what extent this 
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point of view is justified and we consider our result as indication of the impossibility to 
extract the rehable value of from Rr^hadrons ■ 

Authors are thankful to V. L. Chernyak, M. E. Pospelov, A. I. Vainshtein and 
A. S. Yelkhovsky for valuable discussions. This work was supported by the Russian 
Foundation for Fundamental Research under Grant 95-02-04607a. The work of S.F. has 
been supported by the INTAS Grant 93-2492 within the program of ICFPM of support 
for young scientists. 



6 Appendix. Integration over the SU(3) group 

The average values of expressions polynomial in the elements of the SU{N) matrix 
may be easily found by, for example, graphical method described in |^. However, for 
calculation of the integrals including 6'-functions like / \TrO\^6{h) (|33| ) one has to specify 
in explicit form the parametrization of the matrix and the measure of integration. We 
use the following parametrization of the SU{3) group element: 

/ -P^P:S,S2 - PiP3P^CiC2C3 P^P:CiS2 - P2P3P5S1C2C3 P3C2S3 \ 
U=\ P^P^SiC2~PlP^P^CiS2C3 -PIP^CiC2-P2PaP^SiS2C3 P4S2S3 \, (59) 

V -P1C1S3 P2S1S3 P5C3 / 

where: Ci,2 = cos , si,2 = sin , C3 = cos-ip3 and S3 = sinips, Pn = e*"^". The 
parameters ipm and ipn vary in the range: < ^'1,2 <7r, < ip^ < 7i/2 and < v^n < 27r. 
It may be shown that the matrix (|59D is unitary and unimodular. The measure of 
integration over the matrix group U is given by: 

dU = N Juail det {M)\^/^ , (60) 
Mi, = Tr{U~\diU)U-\djU)) . 

In our case the parameters {i = 1, ...,8) are the three angles ipm and five angles (pn- 
After rather tedious algebra one obtaines: 

dU = NUdipml^difn sin ipi sin ip2 sin^ ip^, cos ip3 . (61) 

The results of the averaging with this measure of the expressions {\TrO\^) and 
{\TrO\'^Re[detO{Tr0^y]) coincide with those obtained by the graphical method p9[ . 
The integration of the expressions including 6'-functions has been performed numerically. 
The results are given bellow: 

{\TrO\^) = - ■ {\TrO\^e{h)) = 1.368; 
5 

{\TrO\^Re[detO{TrO+f]) = ^ ; {\TrO\^Re[detO{TrO+)^]e{h)) = 0.473; 

{A\TrO\^) = ^ ; {A\TrO\^9{h)) = 0.510 . (62) 
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